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Section 1 


INTRODUCTION 

This report summarizes the results of a program performed for the 
NASA Langley Research Center by the Calspan Corporation to study windshear 
and gust alleviation methods with consideration of the application of the 
NASA total energy/energy rate probe to this particular task. 

The study has several different objectives and related phases. One 
was to determine whether or not a successful control system design could be 
obtained for the TCV 737 aircraft using only elevator or only throttle, and 
then using both control effectors in concert. A second objective was to 
assess the usefulness of a direct measurement of the shear or turbulence as 
compared to a strict feedback regulation approach to gust alleviation. Still 
another objective was to test the concepts of alleviation with respect to 
ground speed or with respect to airspeed, sometimes referred to as 4D or 3D 
regulation. If the control system of the aircraft is designed to maintain 
a constant airspeed through a windshear, the alleviation system is said to 
be a 3D system. If ground speed is regulated through the windshear, the air- 
craft presumably will land at the touchdown point both in space and on time 
and is referred to as 4D regulation. The airspeed/ground speed regulation 
approach can be considered only with respect to horizontal velocity. However, 
regulation of vertical velocity only with respect to ground can be considered 
as a reasonable criteria in a wind gust/shear alleviation control system 
design study. 

Infused in each of the windshear/gust alleviation designs and design 
philosophy is the application of the NASA total energy/energy rate probe de- 
scribed later in the body of the report. It was felt that the probe could be 
quite useful in this application for three reasons: 

1. A constant rate of change of energy during a landing is con- 
sistent with the 3D philosophy of vehicle regulation in a wind- 
shear. Constant kinetic energy is approximately equivalent 
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to constant airspeed at constant air density while constant 
rate of change of altitude is equivalent to constant rate of 
change of potential energy. It was therefore found appropriate 
to consider energy to be an appropriate criteria in terms of 
direct inclusion in a quadratic performance index criteria for 
design purposes. 

2. The probe as a major component in a control law would be an 
attractive instrumentation feature because of the simplicity 
and versatility of the device. Pneumatic methods of obtaining 
energy rate are straightforward and reliable. Therefore, 
effort was made to include the total energy/energy rate probe 
as a sensor to be incorporated into the control system design 
to simplify the overall system mechanization. 

3. The total energy probe is basically an air data rather than 
an inertia measurement device. Because of this, it was felt 
that the sensor could be useful in the development of instru- 
mentation that would directly detect and measure wind shear. 

The study program was considered exploratory in the sense that sam- 
ple and example methods of the approach to gust and shear alleviation were 
tried rather than trying to take into account all the possible nuances of the 
vehicle or the environment. For instance, secondary effects of a wind shear 
which often occur during a thunderstorm, such as changes in static pressure, 
were ignored as were higher order dynamics and nonlinearities of the airframe. 
The Dryden turbulence model was chosen, partly because of ease of programming, 
but also with the knowledge that other turbulence models, such as the Von Karman 
or Tomlinson models, would have relatively little effect on the results. It is 
clear, however, that a more complete simulation is advisable before flight test 
of any of the several systems described in this report. 

The report is divided into four major sections and three appendices. 
The second section describes the dynamic and kinematic equations applicable to 
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the landing approach flight of an aircraft in wind shear and turbulence. The 
third section describes the control system design approaches used for wind 
shear/gust alleviation, while the fourth section includes the simulation re- 
sults and an analysis of these results. Section five summarizes the results of 
the study and makes recommendations for further investigation of the use of 
the total energy probe in a wind shear/turbulence environment. 

Because the particular computations apply only to the particular 
vehicle used in the study program, the vehicle specific results appear in 
several appendices. These appendices compile the results of the calculations 
performed during this program. 
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Section 2 

DYNAMICS AND KINEMATICS 

2.0 PROBLEM DEFINITION 

In this section, the TCV airplane equations of motion, the shear and 
turbulence models, and the total energy probe measurement model are defined. 

2.1 EQUATIONS OF MOTION 


In defining the equations of motion of the TCV airplane, the follow- 
ing notation is used to identify velocities and angles relative to the ground 
or the atmosphere : 

PA - plane relative to the air. 

PG - plane relative to the ground. 

AG - air relative to the ground. 

It is assumed, in the definition of the equations of motion, that the earth is 
flat and the aircraft is rigid. The nonlinear longitudinal equations of mo- 
tion are defined by the following equations: 


Force Equations - Wind Tunnel Stability Axes 


' v pg - - C D - e ein T PC + cos a PG + T Z. sin “pc 1 -i 

^=2 % 


'PG ‘ C L + acs ' ’PG + 1 + Jj [T z. aOS a PG - T x. 3in “PC 1 wf a 


acs y + q + lj T s. cos a pc - T x. Sin “PC 1 


57.3 


1 PG 9 " a PG 


T = T. qos a m qos B m 
x . T . T . 


T = -T . szn a qos g 

Z • V J • J • 


Thrust components in body axes. 
i- each engine 
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Moment Equation - Body Axis 


<? 


§? . - 3 <lSc c + l [z T 

I m I 1 T . x . 

yy i=l yy x x 


x m T ] 
T . z . J 


6 * q 


Zp j distances from C.G. to a reference point for thrust vector 
i i in body axes'- ft. 


h PC = V PC 3in fpG 


( 2 - 1 ) 


c L m ZT 

qS 

c - 5 - 

D qs 

q = I p V PA 


Lift coefficient. Force perpendicular to velocity vector 
in the plane of symmetry. 

Drag coefficient. Force parallel to velocity vector com- 
ponent in the plane of symmetry. 

Dynamic pressure lb/ft 2 
Density of atmosphere slug/ft 3 


V PA 


True airspeed ft/sec 


q = p = r = 0 


The aerodynamic coefficients are functions of angle of attack, flap 
position, etc . . 


c l - f r V V v 

C D = f ( " ; (2-2) 

C = f ( " ) 

m J 

The inertial positions and velocities are computed from the velocities 
with respect to the body fixed inertial frame as follows: 
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Touchdown 

Point 


X T 


Position and Velocity 


V XE ' V PG ° OS y PG 


PG 


V ZE + V PG SZH y PG 


X„ = 


X o + V XE dt X o ~ V PG ° OS y PG 


Z„ = 


Z + dt = Z + 

o ZE o 


r* 


V PG Sin y PG dt 


Inertial Parameters 
; t 

V = V + 

PG o 


V PG dt 


“pff = °0 + i a pc dt 


9 = 9 + 

o 


q dt 


y PG 9 ' a PG 


(2-3) 


where V , o , and 9 are trim values . 
o 3 o o 
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The guidance information is computed as follows: 



* I X E * Z Y 

— ]J X*, + Z*' sin Y q 


Z e " Z E ~ Z R 


Z E “ X E t<m y 0 


Y i ~ 2 Zff 

Y = tori — 

*E 


where: is desired flight i>ath angle 

altitude error 


Y - y. 


* "I f£ y 

um x ~ y 0 
E 


ILS Glide Slope error 


(2-4) 


The velocity and angular motions with respect to the air mass are 
defined as follows: 
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and Y Pd equations expressed in terms of V^ n and y 


'PA 


PG 


, _ „ -i V P0 Uin 'PG - h AG 

: PA I'yg euc t u 


PG /1G 


PG 


a PA 6 “ Y P/1 


V PA = ^PG COS Y PG * U /1G^ 2 + ^ 7 PG St,?1 Yd ^ “ ^/i^ 2 


PG 4G J 


where : u *» + ve headwind 

AG 


h. _ •v* + ue updraft 

Air 


h PA h PG ~ h AG 


(2-5) 


where and h ^ are the horizontal and vertical wind velocities, respec- 
tively. 


The nonlinear equations (2-1) are linearized about a nominal tra- 
jectory representing the steady state flight along a 3° glide slope. The 
linearized small perturbation equations of motion are expressed in the body 
axis system. The equations of motion are defined as follows: 


A u = X A u + X Aw + (X - W )q - g cos 0 A0 + X. 6, + X s 6 + X £ 6 + X u 

PG u PG w PG q o ^ y o S , t 6 e 8 sp u W 

^ t e sp r 


+ X w 
w w 


to PC ■ V“fc + V“pc + <Z q + V 1 + * sin V s0 + hjt + Z S S e + Z 6 5 sp + Z u u w 

n t e sp 


-h Z W 
w w 


= M U LU PG 


M w LW PG 


M q 
<7 


M a A0 
0 


+ 6 , f IK, 5 + M, 


sp 


sp 


M u 
u w 


+ M W 
w w 


e 


- q 


( 2 - 6 ) 
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where 


A U 
AW 


PG 

PG 


A0 


“Pff - U o 
W PG - W o 

e - e 


C2-7) 


U Q and f/ are trim velocities (air or inertial) along the X and Z body axis, 
u Pq and w pq are the total velocities along X and Z body axis, 9^ is the trim 
pitch attitude, u and w ;) are the wind velocities along the X and Z body axis. 


w 


6 , 6 , and 6 are the throttle, elevator and spoiler control inputs. The 

o & S jp 

perturbation in velocities with respect to the airmass along the X and Z axes 
is defined by 


A u 
Aw 


PA 

PA 


kU PG + U u 
ku PC + 


(2-8) 


where u and w are the wind velocities along the X and Z body axes, respectively. 
W W 

The velocities along the X and Z body axes are defined in terms of the velocity 
in stability axis system by the following equations: 


u 

w 


PG 

PG 


= V. 


PG 


= V. 


PG 


oos a pG 
sin a. pG 


V a 
PG PG 


(2-9) 


where a-^ is the angle of attack. The stability and control derivatives of 

rCr 

Equation (2-6) are given in Appendix A. The equations of motion can be written 
compactly in state variable form as 


x = Fx + Gu + J w 


x = 


A 


- 

PG 


t 

AW„_ 

. u = 

5 

PG 

j 

e 

q 


6 sp 

A0 


* — J 


u 


w 


u = 


w 


w 


J . . — p • • 

'Z'J 


( 2 - 10 ) 
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i - l, 4 

3 = h i 2 


2.2 


SHEAR AND TURBULENCE MODELS 


The wind shear models used in this study are the severe Kennedy inci- 
dent (TOC), (Ref.l ), Philadelphia (T-25A) and the moderate Tower (T-9A) pro- 
files. These profiles define the horizontal and vertical wind velocities, 
and respectively, as a function of altitude. The sign convention is 
defined as: 

u AQ = + ve for headwind (i.e. - X direction in aircraft axes) 

Horizontal component of air motion relative to the ground. 

h.„ = + ve for updraft (i.e. - Z direction in aircraft axes) 

AG 

Vertical component of air motion relative to the ground. 

It is recommended in Ref.l that the shear be defined both as a 
function of altitude and distance. However, in this study the shear profiles 
used are only a function of the altitude. 


The horizontal and vertical velocities are defined by the Dryden 
spectra (Ref. 2 ). The Dryden spectra are defined as: 


C-! - 

<p (Q) = o z * — - — 77 - ; 1 

u u tt 1 + ( L fi) ' 

W u 


where 51 = — 


L 1 + 3 (L 51) 2 

. / n | 2 w a? 

~ °W tt JTTTlIDY 


where is the trim air or groundspeed. The turbulence velocities are obtained by 
filtering white noise using the following transfer functions obtained from 
the Dryden spectra. 
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T (a) 
u 

w 



1 + -y- s 

v 0 


T, (s) 
w 

w 



( 2 - 12 ) 


where the turbulence intensities, and a y and the other parameters , are de- 
fined as follows: 


o =5.1 ft/sec Moderate 
w 

= 7.6 ft/ sec Severe 

a 

w 

a S — 

U (0.177 + 0.000823 

a =o 
u w 

- 1000 
u 

H f t 

£ u' 10 

£ » 1000 ft 

U 

L = 

■' (u.r/y + a.000023 h) 1 - 2 • 


h < 1000 ft 
h > 1000 ft 
h > 1000 ft 
10 < h < 1000 ft 
h < 10 ft 

(2-13) 

h > 1000 ft 
10 < h < 1000 ft 


L = 100 ft h < 10 ft 

u 

Experience at Calspan with the use of Dryden spectra for in-flight 
simulation has indicated the necessity of filtering the turbulence signal. 

This is done to remove the fairly significant low frequency content of the 
signal. In this study, the shear accounts for the low frequency wind distur- 
bance. The shear profiles and the time domain equations to generate the turbu- 
lence are given in Appendix A. 
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2.3 


ENERGY PROBE 


A simple energy probe was developed recently (Ref. 3) by 0 . W. Nicks. 
The energy probe gives a measure of the total energy with respect to the air 
mass. It has been flight tested (Ref. 4) to verify its measurement and recent 
studies (Ref. 5 ) have considered the incorporation of energy rate feedback into 
the control law for control of aircraft in a wind shear environment. 


The energy sensor measures static minus dynamic pressure. If this 
is differentiated, it gives a rate of change of pressure which is related to 
the rate of change in total energy relative to the air mass. The following 
equations model what the probe actually measures: 


' probe 


= P 


statzc 


- q 


P s (h) - ± Pi* - P g (h) 


— 0 V 2 
2 p o IAS 


dP , dp , 9 dV m , 

probe _ s_ ah_ _ 2 „ T 1 


dt 


dh dt 2 pV T. dt 2 v T dh dt 


- V„ 


dp dh 


dP 


or = 


s dh 


dV 


*» T T 

dh dt p o IAS dt 


IAS 


(2-14) 


where is the true airspeed (ft/sec), V is the indicated airspeed (ft/sec), 
p is the density of atmosphere (slug/ft 3 ), is the density of atmosphere at 
sea level and h is the altitude. The pressure measurement is differentiated 
pneumatically using a one second time constant and a transducer converts the 
rate of change of pressure to an electrical signal proportional to energy rate. 


The airspeed can be calculated from the inertial speed and the assumed 
wind speed. The variation of the static pressure with altitude (below 4,000 ft) 
is approximated by the equation 

P = 2115 - ,07196h (2-15) 

s 

The air density and the temperature vary in a thunderstorm as a func- 
tion of position. These variations are given by the following equations: 
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IAS 


P 

p IS' 


Vf 

.000-1383 -rjr .+ l 


3 ‘- if 


1/2 


p = -6.519xl0~ 8 h + .0023763 - Q.OxlO~ 6 AT 


where AT = T -15 
c 


T <v, Temp. in Centigrade 
h 4 4000 ft 


T. = 001975 h + 288.15 + LT Deg. Kelvin 

h 4 4000 ft 


ft/see 

slugs/ft 3 


(2-16) 


In Section 3, a control system based on total energy considerations 
will be described. The energy sensor which has been verified for measurement 
of energy rate can be used in the control system by deriving the energy infor- 
mation from the rate measurement. 

2.4 ESTIMATION OF WIND VELOCITIES 

A technique for estimating the horizontal and vertical components of 
the wind velocity is suggested in this subsection. These velocity components 
can be estimated from the airspeed and the inertial speed. The earth and body 
axis system and the various angles are shown in Figure 1. The notation PG 
represents plane to ground, PA represents plane to air, and AG represents air 
to ground; X ^ , 2^ refer to the body axis system and X & , refer to the earth 
fixed axis system. 



Figure 1. EARTH AND BODY AXES SYSTEMS 
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The horizontal and vertical wind velocities are given by the equations 


V AG - [ V PA aos '‘PA + 

Headwind + ve 

i - - V PG C0s Y PC 


h AG = -[ V PA Sin Y PA - K] 

Updraft + ve 


Z 

e 


- r PC Hn 


Y 


PC 


r PA 6 b - a PA 

(2-17) 

a-. = a v J 

TA vane 

eg 

• • 

where X and Z are horizontal and vertical components of the inertial velocity 

The true airspeed is obtained from the airdata system, the angle of 
attack from the vane corrected for position and flow distortion. The inertial 
velocities are estimated from the DME, MLS, attitude and accelerometer measure- 
ments . 


The accelerations along the X and Z body axes are given by 


x h =n - g sin 0 

D A -7 

(2-18) 

z, = n n + q cos 0 
b Z b 


where n % and are the accelerometer readings in ft/sec. The accelerations 
along the earth X and Z axes are computed using the Euler angle transformation 


— — 


- 



X 

e 


aos 0 sin 0 


n - g sin 0 
% 

z 

e 


-sin 0 aos 0 


n 7 + g aos 0 



— 

~ a b 


(2-19) 
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The inertial position and velocity are estimated by complementary 
filtering of MLS position and aircraft accelerometer signals. The equations 
relating x and z to radar range and elevation are 

x — R cos 0 r (2-20) 

MLS 

z = R sin e MLS 


Smooth estimates of earth- referenced position and velocity are ob- 
tained through blending of the MLS derived data with inertial data using comple- 
mentary filters. The inertial data is obtained from body axis accelerometer 
signals transformed into earth-referenced accelerations as given by Equation 
(2-19). The block diagram of the MLS complementary filters is shown in Figure 2. 




r dme cos q mls 
z " r dme sin Q MLS 



X 


where x 3 z — MLS derived position 

x 3 z = Earth-referenced acceleration 
x 3 z 3 x 3 z = Filtered estimates of earth- 

referenced position and velocit' 
respectively 


Figure 2. MLS COMPLEMENTARY FILTER BLOCK DIAGRAM 
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The equations for the complementary filters are given as follows: 

? (s-hkj)x + kgSx 

x = — — — 

s z +k^s+k^ 

x + (sk +k )x 

/\ J. u 

X - — 

s z +k^s+k^ 


( 2 - 21 ) 


2r , oj 

n 

w 2 
n 

damping ratio 

undamped natural frequency 
MLS derived position information 
earth-referenced acceleration 
filtered estimate of position 
filtered estimate of velocity 

The same filters are used for the estimation of Z position and velocity. 
These complementary filters were used successfully by Calspan (Reference 6) in 
the X-22A flight program. The damping ratio was selected as 0.7 and the 
natural frequency selected as the ratio of accelerometer noise to position 
noise. It is emphasized that the discussion of this subsection represents a 
suggested approach for estimating the wind velocity components and that the 
complementary filters were not used in the control system simulation. 


where . k ^ = 
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Section 3 

CONTROL SYSTEM DESIGN 


3.1 GENERAL DESIGN APPROACH 

The objective of the control system design is to obtain the most 
effective shear/gust alleviation system based upon the three tradeoffs listed 
below . 

1. Number of controllers - i.e., throttle or elevator alone, throttle 
and elevator, and finally throttle, elevator and spoiler. 

2. Shear and gust alleviation based upon assumptions associated 
with the ability to sense the shear and turbulence directly or 
not . 

3. Regulation or alleviation with respect to the airspeed or ground- 
speed. 

Other considerations are associated with the axis system in which the 
aircraft is described. For instance, regulation could be with respect to y and 
VpA (Vp(j) or to ^ atu * u ' Either is acceptable for one represents the flight 
requirements in a polar axis system while the other is rectilinear axis system 
oriented. 


There are at least three different fundamental approaches that can be 
taken in the design of a gust/shear alleviation system. The most basic approach 
and the one that is most often taken is the one of fundamental regulation of the 
plant; i.e., if the feedback is "large enough," the airplane closed-loop dynamics 
will be of high frequency and statically rigid. It will be difficult to perturb 
the aircraft from its initial flight condition. This approach involves feedback 
in the classical negative sense: any feedback signal is used to attenuate the 

perturbation motion that caused the signal. 
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The second approach is to design a system that is less insensitive 
to turbulence or external disturbances. This is done by effectively destabil- 
izing the airplane to a gust or disturbance input. However, if feedback is 

used to effectively reduce C Q , the maneuverability of the aircraft will be 

a 

adversely affected as well. This approach can, therefore, be useful only when 
the aircraft is in an autopilot mode of operation because this design approach 
cannot differentiate between a commanded input and a gust input. 

The two general methods of gust and shear alleviation described 
above involve feedback around the aircraft and, in fact, the airplane must 
respond to the gusts before the motions are sensed and fed back to the vehicle 
controllers. Therefore, the response of the vehicle to turbulence can be 
reduced, but never totally eliminated. 

The third and most direct way to accomplish the gust alleviation 
task is to directly sense the turbulence component of the air mass relative 
to the aircraft and use this signal to directly drive the control surfaces of 
the airplane in such a way that the controllers produce forces and moments 
that directly counter the forces and moments generated on the vehicle by the 
gusts . 


Design methods used in this study include the first and third methods 
described above. 

3.2 SINGLE CONTROLLER DESIGN 

The single controller system makes use only of the elevator for 
flight path control. The regulator control law is obtained by minimizing 
the performance index 

00 

v = j j (q 1 ^ u \ G + q 2 hhl + r6 l )dt f3_1 - ) 

0 

where is the deviation in the inertial speed and A h is the deviation in 

altitude . 
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subject to the aircraft dynamics 
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and U , W and 0 are trim values. The derivation of A h equation is given in 

O O O -I => 

Appendix A. 


The minimization of the performance index through the use of the Euler- 
Lagrange equation results in the control law 
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(3-3) 


where P is the solution of the Riccati equation and A h is the error in altitude 
given by 

A h = h - h (3-4) 

c v 

and h is computed from the range as 

K c = k h R = I sin (3_5) 

and is the glide path angle. The performance index is minimized for several 
sets of values of , q ^ and r and one set is chosen that produces good regu- 
lation in the sense that the perturbation motions are minimized relative to the 
required control activity. 
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As noted earlier, feedback regulation alone has limited ability to pro- 
duce gust and shear alleviation because the motions of the vehicle due to the 
gusts and shears must be sensed before the force and moment producing devices on 
the vehicle can be activated to reduce the effects of gusts and wind shear on 
the airplane. 


Feedforward control using direct gust and shear sensing is sequentially 
added to the control system by calculating the control deflections required to 
totally counter the forces and moments produced on the vehicle by the gusts. The 
linearized equations of motion are defined as 


where u g = [u y w y ] T 


x = Fx + Gu + Ju 

9 


(3-6) 


The total excitation to the aircraft is zero if the controls are 
activated according to the control law 


Gu + Ju = 0 
9 


u 



J u. 


(3-7) 


For the application described in this report, the control effective- 
ness matrix G is singular, so a generalized solution to the control law of 
Equation (3-7) is required. This can be obtained by considering the minimum of 
the mathematical norm 

I = (uC? + u T /) q (Gu + Ju ) (3-8) 

I q gig 

which yields a control law 


Gu ~b Ju 

9 


u = -(G F qG)~ 1 G T J u g (3-9) 

A more detailed discussion of gust and shear alleviation using the 
method of direct measurements of the gusts and wind shears is given in Refer- 
ences 7, 8 and 9. For several of the applications described in this report, only 
the throttle was driven by the sensed wind shear. 
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3.3 


MULTICONTROLLER REGULATOR DESIGN 


The quadratic performance index of linear optimal control serves as 
a criteria for the regulator design investigated in this program. In all cases, 
it was assumed that the inertial states could be measured directly without 
the use of Kalman filtering or other forms of state estimation. For many 
states, such as pitch angle and pitch rate, this approach appears justified 
because the measurements are made directly with respect to an inertial axis 
system. Air data measurements (such as an angle of attack measurement device) 
sense a combination of inertial (steady air mass) quantities and gusts, or 
environmental quantities. They are then used in combination with the inertial 
measurements to obtain direct estimates of the environmental quantities. 

During this study it was assumed that the aircraft and disturbance 
states are measured with sensors that are noise free. In practice, the assump- 
tion of noise free inertial quantity measurements is a good assumption. States 
derived from air data measurements usually require filtering. The required 
filtering is a function of the individual application and can be done really 
successfully only by experimentation. In practice, it is found that accurate 
knowledge of the vehicle stability and control derivatives is of utmost impor- 
tance. The plant noise is usually attributable to unknown derivatives or 
higher order modes excluded from the original vehicle model equations of motion. 
Filtering should be designed to attenuate the higher frequency noise or "obser- 
vation spillover" excluded from the plant model description. As shown in Sec- 
tion 3.5, robust deterministic observers, based upon noise free measurements, 
can be designed as low pass filters that would have the desired result of atten- 
uating the higher frequency plant noise associated with reduced order modeling. 

The quadratic performance index of linear optimal control represents 
an input-output approach to control system design. In general, the closed-loop 
eigenvalues are a direct function of the transmission zeros of the system and 
the weighting matrices in the performance index. The closed-loop eigenvalues 
asymptotically tend toward these transmission zeros as the output weighting 
becomes large with respect to the input (Reference 10) . The response to a 
command input would then tend to resemble the response of a Butterworth filter. 
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In the gust alleviation application, the input is the turbulence or 
wind shear vector, not the actual controllers. Therefore, the optimal response 
of the system should be Butterworth in behavior with respect to a turbulence 
input rather than a command. The problem then reduces to one of defining a per- 
formance index or optimal closed- loop system that reflects the proper inputs. 

The system is defined by the linearized equations of motion 


x(t) = Fx(t) + Gu(t) -h Ju (t) 

9 

y(t) = Hx( t) (3-10) 


where x^(t) = [A u^a)^ A w^Jt), q(t ), A8(£J], u~ = [6 Jt) 3 6 (t ) ] and u (t) 

Jrijr r(jr is & y 

represents the shear/gust vector. The matrices F, G and J represent, respec- 
tively, the system matrix of dimensional stability derivatives, the control 
effectiveness matrix of dimensionalized control derivatives and the gust effec- 
tiveness matrix. It is assumed that the airspeed is closely approximated by 


A u 
AW 
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PA 


ku PG + U w 
“PC + W w 


(3-11) 


The state Au^ then represents the inertial velocity or airspeed be- 

F(s 

fore the vehicle encounters the gust or shear. The gust or shear represents 
a perturbation input to the system. 


A performance index is defined as 


00 

r 


V = rmn 
u 

9 


l yTqy + » g ru g ]dt 


(3-12) 


that yields a closed-loop system matrix 


x = (F-JK)x + Gu + J w 

o g 


(3-13) 


This closed-loop matrix would define a system whose response is op- 
timal, i.e. Butterworth with respect to a gust or shear input. 
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The system realization problem is to define a feedback control law 
u =-KjX such that the closed-loop system matrix is the same as the system 
matrix defined by Equation (3-13), i.e. find such that 


F-JK = F-GK 2 

K 2 = G~ 2 JK = (G F G)~ 1 G T JK (3-14) 

An alternate method is to define an output matrix H 2 such that 

H 1 (Is-F)~ 1 G = H(Is-F)~ 1 J (3-15) 


thereby guaranteeing that the transfer function matrix of the outputs y 2 with 
respect to the control inputs would be the same as the matrix of transfer 
function y with respect to the gust or shear inputs. The performance index 

oo 

V = min I ( + u ru) dt (3-16) 

u J — ” 

o 

then yields a control law u = -K^x that would produce an optimal response 
with respect to the gust and shear inputs rather than the control inputs, 
although the system would be mechanized in terms of the control inputs. 


3.4 MULTICONTROLLER DESIGN 


It is possible, theoretically, to counter exactly the forces and 
moments produced by the wind disturbance with the three controllers, throttle, 
elevator, and spoilers used in a feedforward control law without feedback. 

The aircraft dynamics are given by the equation 
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where U , W . and 0 are trim values. The feedforward control law is given by 

o o’ o 

the equation: 



(3-18) 


Since with three controllers, independent control over the three 
degrees of freedom of longitudinal motion is possible; the control law of 
Equation (3-18) counters the effects of the wind disturbance exactly. It is 
assumed that the exact measurement of wind velocities is possible. As will be 
shown in Section 4, the spoiler control activity is excessive with the use of 
this control law indicating insufficient spoiler control effectiveness for this 
application. 


The objective of this study is to use the two controllers, the 
throttle and elevator, to alleviate the effects of wind disturbance. 


The two-controller control system is based on the performance index 
that would yield a control law that would produce an optimal response with 
respect to the elevator and throttle control inputs. Four types of control 
systems are designed. These are: 

• Control system 1 consists of a feedback part that regulates the 
inertial speed and altitude and a feedforward part that drives 
the throttle using the sensed gust signals to minimize the devia- 
tion in the inertial speed. 
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• Control system 2 consists of a feedback part that regulates the 
inertial speed, altitude, and altitude rate and a feedforward 
part that drives the throttle using the sensed gust signals to 
minimize the deviation in inertial speed. 

• Control system 3 regulates airspeed and altitude or airspeed 
altitude and altitude rate. 

• Control system 4 is obtained by minimizing the deviation in the 
desired energy profile to fly constant airspeed using the throttle 
and minimizing the deviation in the glide path using the elevator. 

Control systems can be designed either to fly constant inertial speed 
or constant airspeed on the glide slope. The control system structure is shown 
in Figure 3. Figure 3a shows the feedback regulation of either the inertial speed 
or airspeed through the feedback of the appropriate measured quantity. Figure 
3b shows the regulation of the inertial speed through the feedback of the meas- 
ured inertial speed and a feedforward part that uses the measured wind distur- 
bance to drive the control surfaces. The feedback gain K is determined using, 
for example, optimal control method and the gain K j using a generlized inverse 
method. Figure 3c shows the regulation of the airspeed using the measured 
inertial speed and wind disturbance. The feedback and feedforward gains are 
equal and same as the gain X in Figures 3a and 3b. The inertial/ airspeed is 
regulated by the control laws in the presence of both steady and accelerating 
winds . 

Control System 1 


The aircraft dynamics used for the design of the feedback regulator 
are given by: 
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A h = sin Q q u pG - cos 0 o w pG + (U Q cos Q q + W q sin B q )LB 
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A u - INERTIAL OR AIRSPEED 

Figure 3a. INERTIAL OR AIRSPEED REGULATION 



Figure 3b. INERTIAL SPEED REGULATION AND DIRECT 
GUST ALLEVIATION 



Figure 3c. IF INERTIAL SPEED IS MEASURED, 
BUT AIRSPEED IS REGULATED 


Figure 3. INERTIAL AND AIRSPEED REGULATION 
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where U Q , W Q and 0^ are trim values and hh is the error in altitude. 


The feedback regulator is formulated using the control law by mini- 
mizing the performance index: 


00 



(q l hU PG + q 2^ hl + r l' & t + r 2 S e )dt 


(3-20) 


o 

The objectives specified in the performance index are to minimize 
deviations in the inertial velocity and altitude to keep the aircraft on the 
glide path. 


The performance index does not include any wind disturbance terms. 
The control law obtained by the performance index minimization regulates 
against the wind disturbance and minimizes the deviation from the glide path. 


The minimization of the performance index is accomplished through the 
use of the Euler -Lagrange equation. The performance index is minimized for 
several sets of q., r- , i=l,2 and the set that produces good regulation for the 

Lr is 

worst case shear is chosen. The feedback control law is given by 
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where P is the solution of the Riccati equation. The performance index weight- 
ing factors, the gain matrix and the closed-loop transfer functions are given 
in Appendix B. 
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The altitude error, A h, is given by 


Lh = h^-h (3-22) 

where h is the desired altitude and h is the actual altitude. The desired 
altitude is computed as 


h 

o 


= \R 


(3-23) 


where R is the range to touchdown and the gain K ^ is determined from the glide 
slope angle as 


K h = | sin y Q \ (3-24) 

The feedforward control is based on measuring the gust (shear and tur- 
bulance) signals and using them directly to drive the control surfaces to coun- 
ter the forces and moments produced by the gusts. It is possible, as discussed 
previously, to completely counter the forces and moments produced by the gusts 
using the three controllers, the elevator, throttle and the spoilers. 


The feedforward control is formulated using two controllers, the 
elevator and the throttle. There are several ways to formulate the feedforward 
control system. One way is to minimize the cost function 


J = ( Gu+Ju ) ( Gu+Ju ) 
9 9 


(3-25) 


where u is the gust vector, u is the control vector, and J is the gust effec- 
9 

tiveness matrix. The resulting controller is given by 


u = -(G r G)~ 1 G T Ju (3-26) 

9 

A second possible feedforward control law is to use the sensed gust 
signals to drive the throttle and elevator to minimize the deviation in veloci- 
ties along the X and Z axes. A third possible control law is to drive the 
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throttle using the sensed horizontal gusts and the elevator using the sensed 
vertical gust to minimize the deviations in velocities along X and Z axes. 

For the severity of wind disturbance considered in this study, these feed- 
forward control laws would introduce undesirable excursions in angular motions 
of the aircraft. Therefore, the feedforward control system formulated uses the 
sensed gust signals to drive only the throttle. The inertial velocity equa- 
tion with gust term is given as 


Am = X Am„ + X Aw + (X - W )q - g cos 0 + X. J, + I. <5 

PG u PG w PG q o 6 t t <5 c 


+ X u + X w 
u w w w 


(3-27) 


The feedforward control law is given by 
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(3-28) 


where the sensed gust signals drive the throttle directly to alleviate the 
effects due to wind gusts on the inertial velocity along the X axis. The coup- 
ling effects due to other responses are alleviated through feedback regulation. 
The complete control law is given by 
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(3-29) 


The comolete control system is shown in Block diagram form in Figure 4. The 
outer loop determines the nroner controls for minimizing the glide slope error. 
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Figure 4. BLOCK DIAGRAM OF CONTROL SYSTEM 1 
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A “PC ” U PG ~ U o 
hW PG = W PG _ W 
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A6 = 8 - 0 

o 

where U Q , W Q and 8^ are trim values. Derivation of t\h equation is given in 
Appendix A. 


index : 


The feedback control law is formulated by minimizing the performance 
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(3-31) 


o 

The objectives specified in the performance index are to minimize 
the deviations in the inertial velocity, altitude, and altitude rate. Since 
the altitude rate is proportional to the potential energy rate, inclusion of 
this term in the performance index can be considered as minimizing the devia- 
tion in the potential energy rate. The inclusion of the rate term in the per- 
formance index will provide lead information about the deviation from the glide 
path due to wind disturbance. The feedback control law is given by 
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(3-32) 


where P is the solution of the Riccati eauation. The performance index 
weighting factors, the gain matrix and the closed-loop transfer functions are 
given in Appendix B. 
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The altitude error and the altitude rate are given by 


A h = | sin y Q | R, Mi = \ sin y^ | R 


(3-33) 


where E and R are the range and the rate of change of range. The feedforward 
control is the same as the feedforward control of control system 1: 
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(3-34) 


The complete control system is shown in block diagram form in Figure 5 



Control System 3 


The objectives in the design of control systems 1 and 2 were to mini- 
mize the deviation in the inertial speed and the glide path. The performance 
indices used in the design of these control systems can also be used to design 
a control that minimizes the deviation in the airspeed and glide path. The 
two performance indices, the first minimizing the deviations in the airspeed 
and altitude, and the second minimizing the deviations in airspeed, altitude, 
and altitude rate, are given as: 
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(3-35) 


where A u pA is the perturbation in the airspeed. The dynamic equations are as 
defined by Equation (3-19) and Equation (3-30) . The same weighting factors 
used in control systems 1 and 2 are used resulting in the same gain matrices. 

The feedback control law obtained from the first performance index 
is given by 



where A h is as defined by Equation (3-22) . 

It was necessary to change the gain on airspeed to the throttle, the 
first element K ^ of K, as defined by: 

aX aX 

K u - xf Au pa = xf (Lu pg * V u - 37) 

s t s t 

This gain is determined based on the feedforward control law as defined by 
Equation (3-28) with the sign on the gain on u changed. The factor a was 
determined to obtain good airspeed regulation for the severe shear. The feed- 
back control law with the feedback of inertial velocities is given by: 
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(3-38) 
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where the gain on A u to <5 is as defined by Equation (3-37) . The feedforward 
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control law is given by 
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The gains on A and u to throttle are the same. 
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(3-39) 


The feedforward 


control changes the inertial speed along the X-axis proportionately to the hor- 
izontal wind speed and compensates directly for the vertical wind. If the air- 
speed is measured and fed back to the throttle with a sufficiently high gain, 
feedback regulation will be sufficient to alleviate the effects of wind dis- 
turbance. The full control law is given by 
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The final gain matrices and the closed-loop transfer functions are 
given in Appendix B. The control system is shown in block diagram form in 
Figure 6. 


The second performance index in Equation (3-35) includes the altitude 
rate in the performance index. The complete control law including the feed- 
back and feedforward parts is given by 
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Figure 6. BLOCK DIAGRAM FOR CONTROL SYSTEM 3 


The gains on A u and u, to the throttle are the same and as defined 

rir W 

by Equation (3-37) . The performance index weighting factors, the gain matrix 
and the closed-loop transfer functions are given in Appendix B. The error in 
altitude and altitude rate are as defined by Equation (3-33) . The energy 
probe can be used as a measure of h in the feedback loop. 


Control System 4 


Control system 4 is designed in two stages. In the first stage a 
performance is defined with the objective of minimizing the deviation in the 
energy expended by the aircraft with respect to the airmass. The total specific 
energy expended by the aircraft with respect to the airmass is given by the 
expression 

U PA 

E = h (3-42) 


where u^ A is the actual airspeed. 
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The desired energy to be expended by the aircraft with respect to 
the airmass to stay on the glide path is given by 

U o 

E c = h c + 2? C 3 -«) 

where h is the proper altitude to stay on the glide path and U Q is the refer- 
ence airspeed. 

A performance index that minimizes the deviation in the energy using 
the throttle control is defined as 


V = j j ( qAE 2 + r&pdt 


subject to the dynamics 
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The A E equation is expressed in terms of the stability and control 
derivatives. The derivation of the A E equation is given in Appendix A. The 
energy term in the performance index could have been expressed as a function 
of altitude and airspeed. However, this would have resulted in the velocity 
term raised to the fourth power. Including the energy term directly in the 
performance index and defining the dynamic equation for A E will result in a 
straightforward synthesis of the control law with a direct feedback from 
energy. Minimization of the performance index results in the control law 
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(3-46) 


The energy sensor probe measures the total energy rate with respect to 
the airmass as defined by 


E = h + 
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(3-47) 


The energy can be obtained from energy rate by integrating Equation (3-47) with 
the proper initial condition for use in the feedback loop. In the second 
stage, another performance index is defined that minimizes the deviation in air- 
speed and altitude using the elevator subject to the dynamics as defined by 
Equation (3-20) . 
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Minimization of the performance index yields the control law 
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The complete control system is obtained by combining the control laws 
of Equations (3-46) and (3-49) and is given by 
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(3-50) 


where the gain matrices K, and are formed from K ^ and . The feedforward 
control that uses the sensed wind shear and turbulence signals to drive the 
control surfaces is not required. The control system block diagram is given 
in Figure 7. The final gain matrix and the closed-loop transfer functions are 
given in Appendix B. A performance index could have been defined including the 
energy term, the throttle and the elevator to perform an integrated design. 
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However, this would have resulted in undesirable feedback from the energy term 
to the elevator. The throttle can change the energy of the aircraft with re- 
spect to the airmass appropriately by proper change in thrust level. Conse- 
quently, the control system design was carried out separately. 


h 



Figure 7. BLOCK DIAGRAM FOR CONTROL SYSTEM 4 
3.5 OBSERVER METHODS OF SYNTHESIS 

Although much progress has been made, the systems that result from 
the application of modem control theory are generally too complex to be 
implemented practically. Feedback from all the states to each of the control- 
lers is generally required. In an attempt to avoid or sidestep the require- 
ment that all the states be measured, state estimator or observer systems have 
been devised in an attempt to reduce the number of required sensors or to im- 
prove the accuracy of the state measurements. This has led to a new crop of 
problems, including reduced robustness through reduction of gain and phase 
margins as compared to state feedback systems, and an increase in the order of 
the dynamic response of the system which always degrades the flying qualities 
of the airplane. 
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The observer designs are based upon the assumption that the sensors 
are in themselves noisy and accurate measurements of the corresponding motions 
of the aircraft cannot be made. Fictitious plant noise, either white or colored, 
is often added to the system description to improve the system sensitivity. 

In practice, the addition of fictitious plant noise is equivalent in many ways 
to an assumed reduction of sensor noise. But sensor noise in itself is an 
unjustified quantity; sensors used in aircraft have very low noise levels. 

To demonstrate, all one has to do is activate a pitch rate gyro or accelero- 
meter and observe the very low noise level. The noise one usually observes on 
the output of an installed instrument is usually due to higher order modeling 
effects such as structural dynamics and vibration, and this should be represented 
as plant noise. 

Because sensors in themselves are noise free, the observer design used 
in this study is deterministic. These observers depend upon accurate knowledge 
of the stability and control derivatives of the aircraft model, and, in fact, 
the robustness of the observer system depends upon this knowledge. In practice, 
this dependency is not too strong, because the poles of the observers are 
directly related to the zeros of the system transfer function matrix. Knowledge 
of stability derivatives such as accurately obtained in a wind tunnel, are 
often dominant. 

The original observer theory by Luenberger (Reference 11) was developed 
as a state reconstruction technique in which the output of the observer system 
approached, as time increased, the state that is to be reconstructed. The 
theory developed by Luenberger described deterministic observers which are in 
themselves unobservable and do not contribute to an increase in the order of the 
response of the system. In general, the poles of the observer can be chosen 
arbitrarily but this choice usually results in a requirement that the control 
input to the aircraft also is a control input to the observer system. In fact, 
by properly selecting the observer poles the observer system itself can be con- 
siderably simplified, as is shown in this section. 
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Consider the linearized, small perturbation equations of vehicle 

motion 


x(t) = Fx(t) + Gu(t) (3-51) 

y(t) = Hx( t) 

where x is the state vector, u the control input vector and y the measurement 
set such as a pitch rate gyro and accelerometer. The observer is of the form 

z(t) = Az(t) + By ( t) + Bu( t) (3-52) 


subject to the conditions : 


1. The matrix A is chosen such that the observer is stable. 


2. A transformation T can be found such that 

TF-AT = BH (3-53) 

3. TG - D where F G and H are defined above. (3-54) 


In order for the observer to be as simple as possible, not requiring 
measurements of the input to the aircraft but operating entirely on the system 
outputs, the following conditions must also be satisfied: 


4. The observer poles, defined from | Js-A | = 0 are chosen from the 

transmission zeros of the system, defined by | H(Is-F) ^G j = 0 

(3-55) 

5. TG = D = 0 (3-56) 
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6 , 


The matrices H and T constitute a transformation M 



Mx 


(3-57) 


such that the equivalent output and observer control law is 
defined by 


u = -Kx = -KM 


rl 



(3-58) 


The condition that TG = D = 0 is automatically satisfied if the ob- 
server poles are chosen from among the system transmission zeros. Not only 
does this simplify the observer network but also reduces the number of required 
measurements of the vehicle dynamics. Condition 6 implies that the output 
measurements must be such that the aircraft dynamics are completely observable 
through the output. This condition is usually easily satisfied with an accel- 
erometer or angle of attack sensor. 


Example : 


Observers can often be used to enhance the observability of the sys- 
tem in the sense that much lower feedback gains are often required to obtain 
the same closed- loop dynamic behavior as compared to state feedback. This comes 
about because the feedback through the observers is shaped as a function of 
frequency and this can often be used to the advantage of the designer. Consider 
the system 



'*1 " 


~o 2 -r 

X 1 


“ 1 ' 



Xr. 

2 

= 

1 6 -12 

X 2 

+ 

0 

u 




_1 6 -IX. 



_ o- 


with a single measurement y = Xj - 

2x ^ 


(3-59) 
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If it is assumed that it is for some reason desirable to alter the 
open-loop characteristic polynomial D(s) = |js-Fj = s 3 + 5s 2 + 5s + 2 to a 
closed-loop characteristic polynomial A (s) = \Is-F+GK\ = s 3 + 5s 2 + 6s + 7 } 
then the state feedback control law is given by 



6x Z 


(3-60) 


According to condition (4), the poles of the observers can be chosen 
from among the system transmission zeros which, for this example, are the 
zeros of the transfer function 

y/u(s) = sl+6s+8 = (_s+2)_ (8+4) (3-61) 

s : ‘+5s 2 +5s+2 s 2 +5s 2 +5s+2 

This output transfer function represents a minimum phase, completely 
observable measurement so two output observers can be constructed with observer 
poles equal to the zeros of the output transfer function; i.e. s^ = -2 3 s^ = -4. 
The observers can then be chosen automatically as 


s 2 = -4z 1 + y 
Z 2 = ~ 2z l + y 


(3-62) 


The matrix T of condition 6 listed above and the control law is most 
easily obtained by first transforming into the phase variable form 


w = F w + G u 
o o 

x = Nw 


(3-63) 


where 


F o = 


~0 

1 

0 “ 


~ 0 ~ 


” 6 

5 

1 ~ 

0 

0 

1 

II 

0 

N = 

-1 

1 

0 

-2 

-5 

- 5 _ 


- 1 - 


- 0 

1 

0 - 


(3-64) 


The transformation = Mx can then be determined as follows (Reference 12): 
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a) The elements of the first row of the matrix M are made up of the 
coefficients of the numerator polynomial of the output transfer 
function, i.e. 


m 1 = [8 6 1} 


(3-65) 


b) The second row of the matrix is made up from the coefficients 
that result when s+4 has been factored out of the numerator 
polynomial of the output transfer function, i.e. 


m^( s ) 


s 2 +6s+8 

s+4 


s+2 


(3-66) 


so 


m 2 = [2 1 0\ 


c) Similarly, the third row is formed from the coefficients that 

result when the factor s+2 has been factored out of the numerator 
polynomial 


m^(s) 


s 2 +6s+8 

s+2 


s+4 


so 


m- = [4 1 0 ] 

fJ 


(3-67) 


tion 


This process is exactly the same as defining the matrix T in the equa- 


TF - AT = BH 

when F and H are defined in the phase variable form. 

The observer control law is then given by 

= J Z 1 - I Z 2 ( 3 - 68 ) 


u = -KM 


-1 


y 

z 
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The closed loop system defined by the control law of Equation (3-68) has the 
same robustness or gain and phase margins as the state feedback system, but the 
system employing observers uses only one output or sensor, the feedback gains 
are much lower and the observers constitute low pass filters, useful for struc- 
tural mode stabilization or noise suppression. 

It is shown above that if the observer poles are chosen from among 
the zeros of the transfer functions, no plant input measurements are required 
in the observer synthesis for systems that are completely observable. For 
multicontroller systems, exactly the same principle holds except that the ob- 
server poles are chosen from the transmission zeros of the system, defined as 
the roots of the polynomial 


\H(Is-F)~ 1 G\ = 0 (3-69) 

when the system is defined by the standard equations of motion 

x(t) = Fx(t) -f Gu(t) 
y(t) = Hx(t) 

As indicated by Equation (3-69), the transmission zeros are a func- 
tion of the sensors or outputs as well as the controller inputs. In general, 
a robust output observer can be obtained if the system is completely controll- 
able and observable with the chosen sensors and if the number of non-minimum 
phase transmission zeros is equal to n-m, where m is the number of sensors used. 
For an nth order system with m independent output measurements and p inde- 
pendent inputs, the maximum number of transmission zeros is given by 

* = ^ (m-pTl p! C 3 - 7 °) 

Ordinarily, the number of transmission zeros exceeds the number of observers 
to be constructed so there is design freedom, as shown by the following example. 

Consider the 4th order system given by 
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(3-71) 


According to Equation (3- 70) , the maximum number of transmission zeros of this 
system is 


N = ( n-p) 


r- — 


- — 

m! 

= (4-2) 

3! 

_ (m-p) ! p! 

L 11 21 J 


(3-72) 


A maximum of two of a possible six transmission zeros must be minimum 
phase in order to be able to construct a stable output observer system. There 
are three sets of transmission zeros that can be used in the selection of the 
output observer poles. These three sets involve the outputs (y 2 y p , iy 3 ) 
and y 2 ,y 3 ) and produce the transmission zeros 


(y-L y 2 ) -3s 2 -20s-31 = -3( s+2 . 45) (s+4 .215) 

(y 2 y 3 ) 3s z -8s+6 = 3 ( 3 - j (3-73) 

(y 2 y 7 ) 2s+6 = 2(s+3) 


The system possesses five of a possible six transmission zeros and three 
are non-minimum phase. Any two can be chosen for the system design. The obvious 
choice is to design an output observer using the two measurements y 2 and y 
The third measurement y 7 is not required. In fact, an observer set for the 

tJ 

measurement pair y 2 y 2 is given by 



(3-74) 
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Appendix C defines a multicontroller observer configuration for the 
multivariable control system 1 defined previously in this section. It is 
included to show how the control system would actually be mechanized for flight 
test purposes. 


Section 4 

SIMULATION RESULTS 


This section presents the simulation results of the TCV aircraft 
flying in a wind shear and turbulence environment on a 3° glide slope for 
touchdown. The linearized model of the TCV aircraft, linearized for the land- 
ing flight condition, is used in the simulation. The flight conditions (trim 
conditions) are summarized as follows: 


U (airspeed along X-axis) = 213.92 ft/sec W (airspeed along Z-axis) = 8.63 


a (angle of attack) = 2.31° 
(flight path angle) = -3° 
(elevator) = 2.7° 

Flaps = 40° 


ft/sec) 

6 o (pitch attitude) = -.69° 

5 (throttle) = 9000 lb 

6 (spoilers) = 0° 


Two severe shear profiles (Kennedy and Philadelphia incidents) and 
one moderate shear profile (Tower) are used in the simulation. These profiles 
define the wind velocities as a function of altitude. The turbulence is gen- 
erated from the Dryden spectra. This turbulence is filtered to remove the 
low frequency content. The shear represents the low frequency part of the 
wind disturbance. A a of 7.6 ft/sec is used for the Dryden spectra. The 
shear profiles and the generation of turbulence is given in Appendix A. Non- 
linearities and engine dynamics were not included in this simulation. 

The simulation results presented in this section are summarized as 

follows : 

• Controls fixed simulation in a windshear and windshear and turbu- 
lence environment. 

• Simulation with the elevator, throttle and spoilers driven by the 
sensed shear and turbulence signals. 

• Single control (elevator) system simulation with Kennedy incident 
shear. 
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• Simulation of control system 1 (feedback only) with the Kennedy 
incident horizontal shear only. 

• Simulation of control system 1 with both the horizontal and ver- 
tical Kennedy incident shear, the shear and turbulence, and the 
other two shear profiles (Philadelphia and Tower) . 

• Simulation of control system 2 with all the shear profiles. 

• Simulation of control systems 3 and 4 with the Kennedy incident 
shear. 

For the control-fixed simulation the following responses are shown: 
the shear and turbulence profiles, desired and actual altitude vs. range, per- 
turbation in inertial velocity along the X-axis, altitude rate, the flight 
path angle, inertial angle of attack, pitch rate, pitch attitude, angle of 
attack with respect to the airmass, y with respect to airmass, airspeed and 
time histories of shear and turbulence. For the feedforward control simula- 
tion, in addition to the above time histories, the control activities are shown. 
For the remainder of the simulations, the following responses are shown: the 

wind profile, the desired and actual altitude vs. range, perturbation in iner- 
tial velocity along the X-axis, the altitude rate, the flight path angle, angle 
of attack with respect to the airmass, pitch rate, pitch attitude, the air- 
speed, elevator and throttle control activities. All the responses shown are 
total . 


The first set of simulation results shown are the airplane flying in 
a wind shear environment with the controls fixed. The Kennedy incident severe 
shear profile is used in the simulation. The time histories of the controls- 
fixed simulation are shown in Figure 8. 

For almost the first 50 seconds, the vertical shear is practically 
zero. During this period, the horizontal shear has an initial increasing 
headwind component, then steadies out, and at about 35 seconds has a further 
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Figure 8. OPEN LOOP SIMULATION WITH 
(Cont'd) KENNEDY SEVERE SHEAR 
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Figure 8. OPEN LOOP SIMULATION W I TO 
(Cont’d) KENNEDY SEVERE SHEAR 







increase in headwind component increasing to about 37 kts . at its peak. 

During this period, the angle of attack increases, the flight path becomes 
shallower, and consequently the aircraft is above the glide path. Beyond this 
point, the headwind reverses to an increasing tailwind and there is also a down 
draft. The altitude starts decreasing. At about 75 seconds, the combination 
of the increasing tailwind and down draft forces the flight path angle to 
steepen, decreases airspeed and increases the angle of attack, causing a dra- 
matic drop in altitude. In the simulation, the altitude was not allowed to go 
below 50 ft. which is reflected in the altitude vs. range plot. The controls 
fixed simulation demonstrates the severe effects of wind shear on aircraft 
flight during landing. Simulation time histories with both shear and turbu- 
lence are shown in Figure 9. The results are similar to the shears only case, 
except the responses have more high frequency content. 

Simulation time histories of the feedforward control system using 
elevator, throttle and spoilers are shown in Figure 10, The sensed turbulence 
and shear signals are used directly to drive the control surfaces. The ele- 
vator, throttle and spoilers are used to exactly counter the pitching moment, 

X and Z forces produced by the wind disturbance. The control activities are 
such as to counter the perturbations in inertial responses due to shear and 
turbulence. However, control activities are excessive, especially the spoilers, 
but these can be reduced by simply scaling down the gains to the controllers. 

Simulation time histories of the single control system with the Ken- 
nedy incident shear are shown in Figure 11. The elevator is used to minimize 
the deviation in the glide path. Initially when the headwind is encountered, 
the airplane is pitched down to keep the aircraft on the glide path. When 
the headwind turns to tailwind and the severe downdraft is encountered, the 
elevator alone cannot prevent the deviation from the glide path. The airplane 
loses airspeed and altitude. In the simulation, the altitude was not allowed 
to get below 50 feet. The time histories clearly show the instability of the 
aircraft flight in this windshear environment. The elevator can counter the 
forces and moments of headwind but not the tailwind and the severe downdraft. 
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Figure 9. OPEN LOOP SIMULATION WITH KENNEDY 
(Cont'd) SEVERE SHEAR AND TURBULENCE 
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Figure 10. SIMULATION OF FEEDFORWARD CONTROL THAT 
(Cont'd) DRIVES THE ELEVATOR, THROTTLE, AND 
SPOILERS USING THE SENSED SHEAR AND 
TURBULENCE _ 
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Figure 11. SIMULATION OF THE SINGLE CONTROL SYSTEM WITH THE KENNEDY SHEAR 
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Figure 11. SIMULATION OF THE SINGLE CONTROL SYSTEM WITH THE KENNEDY SHEAR 
(Cont 'd) 


The time histories of the control system 1 (feedback regulation only) 
simulation with the Kennedy incident horizontal shear only is shown in Figure 12. 
Initially when the airplane encounters the headwind, the elevator and throttle 
are activated to reduce the angle of attack with respect to the airmass and the 
attitude and to regulate the inertial speed and altitude rate to stay on the 
glide path. The airspeed increases as the headwind increases. The simulation 
was not carried far enough for the headwind to tailwind. 

Simulation of the control system 1 including the feedforward control 
with the Kennedy incident horizontal and vertical shears is shown in Figure 13. 
Initially, when the headwind is encountered, the throttle and elevator are 
activated to minimize the deviation in the inertial speed by increasing the 
airspeed and reducing the angle of attack and attitude. The glide path is well 
regulated as shown by the altitude and flight path angle time history. As the 
headwind turns to tailwind, the throttle activity is decreased and the elevator 
is activated to increase the angle of attack. When the severe down draft is 
encountered, the throttle activity is increased to minimize the deviation in the 
glide path. The angle of attack decreases to a minimum of -4.5° and increases 
to a maximum of about 6.3°. The elevator is deflected about 12° maximum from 
its trim position. The thrust is increased to a maximum of about 11,700 lbs from 
its trim setting. The variation in the inertial speed is minimized and the 

airspeed variation is proportional to windspeed. 

Simulation of this control system with the shears and turbulence 
is shown in Figure 14. The flight path angle is again very well regulated. 
Examination of the control motions indicates that the throttle activity is high. 
This is because the sensed turbulence is fed to the throttle. In practice, 
this level of throttle activity is undesirable and the throttle dynamics are not 
fast enough to respond to the turbulence. A more desirable means of control 
would be to use sensed shear signals, which are low frequency in content, to 
drive the throttle. Turbulence alleviation can be accomplished through the ele- 
vator or spoilers either by feedback regulation and/or direct gust sensing. 
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Figure 12. SIMULATION OF CONTROL SYSTEM 1 
(FEEDBACK ONLY) WITH KENNEDY 
HORIZONTAL SHEAR ONLY 
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Figure 12. SIMULATION OF CONTROL SYSTEM 1 
(Cont'd) (FEEDBACK ONLY) WITH KENNEDY 
HORIZONTAL SHEAR ONLY 
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Figure 12. SIMULATION OF CONTROL SYSTEM 1 
(Cont'd) (FEEDBACK ONLY) WITH KENNEDY 
HORIZONTAL SHEAR ONLY 
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Figure 13. SIMULATION OF CONTROL SYSTEM 1 (FEEDBACK AND 
FEEDFORWARD) WITH THE KENNEDY SHEAR 
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Figure 13. SIMULATION OF CONTROL SYSTEM 1 (FEEDBACK AND 
Cont'd FEEDFORWARD) WITH THE KENNEDY SHEAR 
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Figure 13. SIMULATION OF CONTROL SYSTEM 1 (FEEDBACK AND 
Cont'd FEEDFORWARD) WITH THE KENNEDY SHEAR 
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Figure 14. SIMULATION OF CONTROL SYSTEM 1 (FEEDBACK AND 
Cont'd FEEDFORWARD) WITH KENNEDY SHEAR AND TURBULENCE 
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Figure 14. SIMULATION OF C0NTR( 
Cont'd FEEDFORWARD) WITH K1 
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Simulation results of the control system with the severe Philadelphia 
T-25A shear and the Tower T-9A moderate shear are shown in Figures 15 and 16, 
respectively. As indicated by the time histories, the flight path angle is 
very well regulated. The performance is similar to the performance with the 
Kennedy shear. 

Simulation results of control system 2 regulating the inertial velo- 
city, altitude and altitude rate, and including the feedforward control for 
the three shear profiles are shown in Figures 17, 18 and 19, respectively. The 
flight path angle for all three cases is very well regulated. 

Simulation time histories of control system 3 are shown in Figures 
20 and 21. Figure 20 shows time histories of the altitude and airspeed regu- 
lation control system simulation and Figure 21 shows the simulation of the 
altitude, altitude rate, and airspeed regulation control system. Initially 
when the headwind is encountered, the throttle and the elevator are activated 
to minimize the deviation in airspeed and glide path. The throttle activity 
is decreased because of the headwind whereas control system 1 increased the 
throttle activity to minimize the deviation in inertial speed. The throttle 
activity is increased when the headwind turns to a tailwind and the severe 
down draft is encountered to minimize the deviation in the airspeed. In con- 
trast, control system 1 decreased the throttle activity after getting past 
the severe down draft. The angle of attack increases to a maximum of about 
5.3° from its trim value. The thrust has increased to a maximum of about 
21,000 lbs from the trim setting and the maximum elevator deflection is about 
6.8°. The flight path angle is well regulated. The inertial speed varies 
proportionately with the windspeed in an opposite sense in contrast to control 
system 1 simulation where the airspeed variation was proportional to the 
windspeed. The angle of attack excursion is smaller than that for control 
system 1. The elevator deflection is about half of the deflection for control 
system 1 whereas the maximum thrust is about two times more than that for 
control system 1. The performance of control system 3 regulating the airspeed, 
altitude and altitude rate is similar to the performance of the airspeed and 
altitude regulation system. 
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Figure 16. SIMULATION OF CONTROL SYSTEM 1 (FEEDBACK) WITH 
THE MODERATE SHEAR 
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Figure 16. SIMULATION OF CONTROL SYSTEM 1 (FEEDBACK) WITH 
Cont'd THE MODERATE SHEAR 
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Figure 18. SIMULATION OF CONTROL SYSTEM 2 (FEEDBACK AND 

FEEDFORWARD) WITH THE PHILADELPHIA SEVERE SHEAR 
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Figure 18. SIMULATION OF CONTROL SYSTEM 2 (FEEDBACK AND 
Cont'd FEEDFORWARD) WITH THE PHILADELPHIA SEVERE SHEAR 
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Figure 18. SIMULATION OF CONTROL SYSTEM 2 (FEEDBACK AND 
Cont'd FEEDFORWARD) WITH THE PHILADELPHIA SEVERE SHEAR 
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Figure 19. SIMULATION OF CONTROL SYSTEM 2 (FEEDBACK) 
WITH THE MODERATE SHEAR 
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Figure 19 . SIMULATION OF CONTROL SYSTEM 2 (FEEDBACK ) 
Cont'd ' WITH THE MODERATE SHEAR 
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Figure 19. SIMULATION OF CONTROL SYSTEM 2 (FEEDBACK) 
Cont’d WITH THE MODERATE SHEAR 
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Figure 19. SIMULATION OF CONTROL SYSTEM 2 (FEEDBACK) 
Cont'd WITH THE MODERATE SHEAR 
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Figure 20. SIMULATION OF CONTROL SYSTEM 3 
(Cont'd) (ALTITUDE REGULATION) WITH 
KENNEDY SEVERE SHEAR 
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Figure 20. SIMULATION OF CONTROL SYSTEM 3 

(Cont'd) (ALTITUDE REGULATION) WITH 

KENNEDY SEVERE SHEAR 
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Figure 20. SIMULATION OF CONTROL SYSTEM 3 
(Cont'd) (ALTITUDE REGULATION) WITH 
KENNEDY SEVERE SHEAR 
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Figure 21. SIMULATION OF CONTROL SYSTEMS 
(Cont'd) (ALTITUDE AND ALTITUDE RATE FEED- 
BACK) WITH THE KENNEDY SEVERE SHEAR 
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Figure 21. SIMULATION OF CONTROL SYSTEMS 
(Cont’d) (ALTITUDE AND ALTITUDE RATE FEED- 
BACK) WITH THE KENNEDY SEVERE 
SHEAR 
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The final simulation results shown are for control system 4. The 
time histories are shown in Figure 22. This control system regulates the 
airspeed through energy considerations. As with the control system 3, when the 
headwind is encountered, the throttle activity is decreased and when the 
headwind turns to tailwind and the down draft is encountered, the throttle 
activity is increased to minimize the deviation in airspeed. The flight path 
angle is well regulated. The angle of attack increases to a maximum of about 3°. 
The thrust has increased to a maximum of 20,000 lbs from the trim setting and 
the maximum elevator deflection is about 5.6° from the trim position. The 
angular motions and the control activity are smaller than that for control 
system 3. 
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Figure 22. SIMULATION OF CONTROL SYSTEM 4 
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Figure 22. SIMULATION OF CONTROL SYSTEM 4 
(Cont ' d) 
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Figure 22. SIMULATION OF CONTROL SYSTEM 4 
(Cont 'd) 
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Section 5 


CONCLUSIONS AND RECOMMENDATIONS 


S.l CONCLUSIONS 

The research performed in this program was done in response to a need 
to develop automatic control methods that would enable an aircraft to safely 
fly through a severe wind shear and turbulence environment. As demonstrated 
in this report, there are many different concepts and criteria associated with 
shear and gust alleviation, and the total list of possible ways to do it are 
far from exhausted. 

The objective of this study program was to demonstrate feasibility 
in a preliminary way. Basic studies were performed such as inertial or ground 
speed regulation compared to airspeed regulation. The purpose was to determine 
the control law differences in the concepts, the control power requirements and 
the basic safety of flight parameters, such as angle of attack excursions of 
the aircraft as it flew through the shear environment. It was shown that regu- 
lation with respect to airspeed or groundspeed is conceptually the same* 
depending upon the variable used in the control law, the aircraft can be regu- 
lated with respect to either with equal effectiveness. This applies, of course, 
only with respect to the horizontal shear. Vertical motions of the aircraft 
must be regulated with respect to ground. 

Basically, then, this program involved the study of criteria defini- 
tion for the gust and wind shear problem. It was not found to be difficult to 
formulate a performance index that would produce a shear alleviation control 
law, either with respect to inertial speed or airspeed. One of the better cri- 
teria involved the use of an expression for total energy in the performance 
index. This produced a control law that was , in many ways, directly related 
to the kind of behavior inherently expected of a gust/shear alleviation system. 

The total energy probe was shown to be a reasonable approach to an 
implementation of a gust/shear alleviation system. Total energy can be expressed 
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as a state variable and, therefore, not only expresses a very good criteria for 
shear alleviation, but the energy probe output also becomes part of the feedback 
control law. 

For the most part, optimal control methods were used to obtain both 
criteria and control laws for the solution of the shear/ gust problem. In general, 
these control laws are quite complex, requiring feedback from all of the states 
to each of the controllers. A theoretical contribution of the research reported 
upon involves the development of robust observers. The application of these 
methods should result in simplified control laws involving requirements for 
fewer sensors. 

5.2 REC OMMEND AT IONS 


1 . The research results reported upon in this report demonstrate the 
feasibility of automatic gust and shear alleviation systems for commercial 
aircraft such as a Boeing 737. In order to more realistically verify the direct 
suitability of the systems defined, it will be necessary to demonstrate their 
usefulness on a simulation that more realistically reproduces the actual envir- 
onment. Among the more important nonlinear and dynamical effects that should 
be added are: 

a) Throttle and control surface servo dynamic behavior. The non- 
linear engine response characteristics are probably the most 
important effects not included in the simulation of the present 
study. 

b) Aerodynamic nonlinearities. 

c) The energy probe mathematical model should be more accurately 
represented. 

d) Sensor dynamics and other characteristics such as biases, mis- 
alignments and noise, if determinable, should be added to the 
system simulation. 
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e) Other turbulence models, such as the Tomlinson representation 
of turbulence, should be tested. 

2. One of the important findings of the present program is the conclu- 
sion that the spoilers of the example aircraft are of limited effectiveness 
and control power. The use of flaps as active devices for gust/shear allevia- 
tion should be explored. 

3. The investigation of suitable performance indices or criteria should 
be continued, the possibilities are far from exhaused. In particular, the 
use of energy rate, E, should be investigated. 

4. Several theoretical areas show potential in producing more effective 
and more realistic gust alleviation systems. The realistic problems of instru- 
mentation should be included in a further investigation of direct gust allevia- 
tion methods. As now formulated, the integral performance criteria only in- 
directly express alleviation criteria, the minimization is with respect to the 
control inputs rather than the environmental disturbances . 

5. More emphasis should be placed on the use of the energy probe in 
terms of criteria definition, control law mechanization and environmental 
disturbance sensing. 

6. A sensitivity analysis of off-nominal flight conditions should be 
conducted to verify the robustness of the control system design. In particu- 
lar, the final control system configuration should be capable of accommodating 
both frontside and backside power-required operation. 
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Appendix A 

EQUATIONS OF MOTION AND WIND MODEL 
EQUATIONS FOR SIMULATION 


This appendix gives the equations of motion of the TCV airplane used 
in control system design and the wind model equations for use in the simulation. 
The linearized aircraft equations of motion used in this study are represented 
in body axis by the following equations: 
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and w w are the wind shear and turbulence along the X and Z body 
axes, respectively with the sign convention that headwind and updraft are pos- 
itive . 


The matrices are defined in terms of the dimensional derivatives as 

follows : 
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The expression for A h is determined as follows: 
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where U-nn an <l w -nn are total inertial velocities along the X and Z body axes, 

r Lr r(j 

respectively; jj W and 6^ are trim values. 
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using small angle approximations, the expression for A h reduces to 
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The expression for Eh is obtained by differentiating Equation (A-5) 
and is given as 
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The rate of change of energy deviation is defined as follows: 
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where u pA is the total airspeed along the X-axis. Since U Q , the trim airspeed, 
is constant 
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The equations of motion (A-l) represent the perturbation in air or inertial 
speed in the absence of wind disturbance. Using the approximation U PA~~ U o ■ 
the AZ? equation can be written as 
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NUMERICAL VALUES OF STABILITY AND CONTROL DERIVATIVES 


Units 

hu pG , kw pG - ft/ sec 
q, A0 - deg/sec, deg 
8 ± - lbs. 

8 g - degrees 
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e = -.69 deg 

o 



Derivative Value 


u 


W 


X 

9 

\ 

h 

X & 
X 6 

z S P 

u 


w 


.0376 

.106 

0 

.00038 

.0065 

-4.37x10 

- .278 

- .711 


A- 3 


Derivative 


Value 


■2.99 E-7 


6 

M S P 


8.12x10 

-.011 

- .361 

- .523 
-3.27 E-4 

3.61 E-4 
-1.215 


Turbulence Simulation Equations 


The transfer functions of the turbulence are given as 
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where w ^ and are two independent white noise sources, 


The time domain equations, using the phase variable representation, 
can be written straightforwardly as 
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The turbulence signals generated from these equations have significant 
low frequency content. In this study the wind shear represents the low frequency 
wind disturbence. Consequently during the simulation of the control system, the 
turbulence signals generated from Eq. (A-ll) were filtered to remove the low 
frequency content using the following filter: 


H(s) = 


,8(s + .25) 

(s + 3.12) (s +10) 


(A-12) 


where the frequencies are in the units of radians. 

The three windshear profiles useddin this study are shown in Figures 
A-l, A- 2 and A-3. 
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Figure A- 2. PHILADELPHIA (T-25A) SEVERE 
SHEAR PROFILE 
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Figure A- 3. MODERATE T-9A SHEAR PROFILE 
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Appendix B 

CONTROL SYSTEM GAINS AND TRANSFER FUNCTIONS 


This appendix gives the weighting factors used in the performance index, 
the gain matrix, and the transfer functions for each control system. 
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Note : The last row of the F and G matrices represents the A E equation. The 

elements of these rows were computed from Eq.(A-8) using the stability 
and control derivatives supplied by NASA. The coefficient of A6 (element 
F(7,4)) in Eq.(A-9) is negative for the particular flight condition con- 
sidered, whereas the coefficient computed from Eq.(8) is positive. This 
discrepancy is because of the approximations involved in linearization 
of the nonlinear aerodynamic equations in obtaining the stability and 
control derivatives, which were used in computing the coefficients of 
Eq.(A-8), and should have little effect on the control law design for 
system 4. 
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Single control system (elevator) 


Performance index to regulate inertial speed and altitude 
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Control System 1 


Performance index is defined to minimize the deviations in the iner- 
tial speed and altitude. 

CO 
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v = ~2 (10bM* G +l00kh 2 + .162 + &*)dt 
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Total control law 
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Control system 2 


Performance index is defined to minimize the deviations in inertial 
speed, altitude and altitude rate. 

CO 

V = 4- f (lOAu 2 ^ lOOAh 2 + 100 Ah 2 + .16 2 + S 2 )dt 

2 rir V Q, 

O 

Total control law 



Units 

A u pG , Aw pQ , Ah - ft/ sec 

q i A0 j S g - deg/sec, deg, deg 
- lbs 
Ah - ft 
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Control system 3 


Performance indices are the same as for control systems 1 and 2. 
The gain on the velocity is changed to 

= ~v r ~~ = -5.03 (99.4) = -500 
6 * 

Total control law to minimize airspeed and altitude deviation: 



Total control law to minimize airspeed, altitude, and altitude rate 

deviation: 



Units 

, tM pG j hh - ft/sec 
q , A9, - deg/sec, deg, deg 

6 t - lbs 
th - ft 
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Control system 4 


Performance index to minimize deviation in energy using the throttle: 



(100EE 2 + 5 2 )dt 


o 

Performance index to regulate altitude and airspeed: 


7 f 

7*4-1 ( lOOObu 2 + lOOhh 2 + S 2 )dt 
2 } PA e 

o 

The gain on A E obtained from the first performance is replaced by a 
higher gain obtained from a higher weighting. The final control law is given 
by 
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Appendix C 
OBSERVER DESIGN 


This appendix describes the application of the observer system design 
methods described in Section 3 of this report. The control law resulting from 
the regulation of Au and Lh using both throttle and elevator is described mech- 
anized as an observer system. Three sensors were assumed to be available for 
feedback control purposes, one for each degree of freedom of motion of the air- 
craft. In general, this is the maximum that would normally be required to in- 
sure system realizability for most physical plants such as aircraft. The three 
sensors chosen were b.u(t) , A h(t) and q(t) . Only the A h(t) measurement is 
absolutely required, others could have been substituted for the A u(t) and q(t) 
states. The altitude error, however, is unobservable in the control sense 
through any other single sensor, so it is a required measurement for the observer 
system. All three sensors represent inertial measurement and are obtainable 
with relatively low inherent sensor noise content. A deterministic observer 
system is considered a realistic way to define a control system for this appli- 
cation to preserve the robustness of the state feedback shear alleviation con- 
trol law. 


Among the groups of measurements tfo and u, Lh and q, and A u and q, 
there are three minimum phase transmission zeros located at s = -4.525, 
s = -2.39 and s = -0.6497. The zeros at s = -4.525 and s = -2.39 were chosen 
for the observer design to replace two state measurements because they would 
represent more stabel observer poles and because this selection resulted in 
observer parameters that could easily be mechanized. Applying Equations of 
Section 3 yielded the observer 
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(C-l) 

When this observer was substituted into the control law of Control 
System 1 to replace the state measurements Aco and A0, the resulting control 
law is given by 


C-l 


A u( t) 


S t ( t) 


.981 

.779 

.025 

- .241 

— 

- .219 

A u( t) 

_V W J 


-2.42 

-.801 

-2 . 966 

14.28 

21.52 

A h( t) 
q( t) 


Zj(t) 

Z 2 ( t) 


(C— 2) 


The observer control law produced somewhat lower overall feedback 
gains, indicating that the observers had not reduced the observability of the 
system and, in fact, had improved it somewhat. It is felt, therefore, that the 
observer system should be easily mechanizable. 
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